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Motion Artifacts in MRI: A Complex 
Problem With Many Partial Solutions 


Maxim Zaitsev, PhD,'* Julian Maclaren, PhD,''* and Michael Herbst, PhD? 


Subject motion during magnetic resonance imaging (MRI) has been problematic since its introduction as a clinical imag- 
ing modality. While sensitivity to particle motion or blood flow can be used to provide useful image contrast, bulk 
motion presents a considerable problem in the majority of clinical applications. It is one of the most frequent sources of 
artifacts. Over 30 years of research have produced numerous methods to mitigate or correct for motion artifacts, but 
no single method can be applied in all imaging situations. Instead, a “toolbox” of methods exists, where each tool is 
suitable for some tasks, but not for others. This article reviews the origins of motion artifacts and presents current miti- 
gation and correction methods. In some imaging situations, the currently available motion correction tools are highly 
effective; in other cases, appropriate tools still need to be developed. It seems likely that this multifaceted approach 
will be what eventually solves the motion sensitivity problem in MRI, rather than a single solution that is effective in all 
situations. This review places a strong emphasis on explaining the physics behind the occurrence of such artifacts, with 


the aim of aiding artifact detection and mitigation in particular clinical situations. 


Ce with other imaging modalities, such as ultra- 
sound or computed tomography, MRI has always been 
particularly sensitive to subject motion. This is primarily 
due to the prolonged time required for most MR imaging 
sequences to collect sufficient data to form an image. This 
is far longer than the timescale of most types of physiologi- 
cal motion, including involuntary movements, cardiac and 
respiratory motion, gastrointestinal peristalsis, vessel pulsa- 
tion, and blood and CSF flow. The effects of motion have 
been well known since the early days of MRI and include 
blurring and ghosting in the image.'” 

Recent technological improvements have improved the 
situation in some cases, but have exacerbated it in others. 
On one hand, incremental performance gains in hardware 
(e.g., higher performance gradients) have enabled faster 
imaging, as have breakthroughs such as parallel imaging. 
Faster imaging means that some scans can be performed in 
a shorter time, leading to a smaller chance of involuntarily 
subject motion. On the other hand, hardware improvements 
have improved the achievable resolution and signal-to-noise 
ratio (SNR), and the sensitivity to motion has, therefore, 
increased. Stronger gradients also mean greater phase accu- 
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mulation for moving spins. Higher main field strengths and 
stronger gradients also mean that a typical MR scan is gen- 
erally louder than it was in the early days of MRI, which 
reduces the chance of infants sleeping through the proce- 
dure. This is particularly true for diffusion-weighted or dif- 
fusion tensor imaging, as was reviewed recently by Le Bihan 
et al.? 

At present, there is no sign that the problem of subject 
motion during MRI examinations will be resolved through 
hardware improvements. The potential of accelerated imag- 
ing seems to be increasingly limited by biologic constraints: 
peripheral nerve stimulation limits gradient switching 
speeds; specific adsorption rate (SAR) limits the use of RF 
excitation pulses; and T1 and T2 relaxation times constrain 
the sequence repetition and echo times, depending on the 
required contrast. 

Furthermore, it seems increasingly likely that there is 
not, and will not be, a single methodological solution to the 
problem of motion in MRI; rather, a toolbox of solutions 
exists. The choice of which tool to apply depends on the 
task at hand, namely, the imaging situation and type of 
motion Some correction 


in question. prevention or 
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techniques are highly effective in particular circumstances, 
but useless in others. Other techniques are relatively general, 
but are not completely effective. It is, therefore, extremely 
important to recognize and understand motion artifacts to 
identify the physical cause of the problem in a particular 
clinical situation and use the best possible tool for the 
problem. 

In this review, we first examine the physical origins of 
motion artifacts and then summarize the main motion pre- 
vention and correction tools currently in use or in develop- 
ment. Our aim is to provide a condensed guide to the 
function, application and limitations of these methods. 


Physical Principles Behind Motion Artifacts 


The appearance of motion artifacts in an image is a result 
of a complicated interaction between the image structure, 
the type of motion, the specifics of the MR pulse sequence 
and the k-space acquisition strategy. In this section we 
briefly revisit the most important physical mechanisms 
affecting the appearance of motion artifacts. 


k-Space and the Image Acquisition Process 

Spatial encoding in MRI is an intrinsically slow and sequen- 
tial process. To appreciate the effect of motion on MR 
images it is important to realize that the data acquisition 
occurs not directly in image space, as is the case in photog- 
raphy, but rather in frequency or Fourier space, which is 
commonly termed “k-space”. 

k-space corresponds to the spectrum of the spatial fre- 
quencies of the imaged object and depending on the imag- 
ing situation can be two- or three-dimensional. Spatial 
frequency spectra of objects with small number of contrast 
features and with smooth intensity variations are predomi- 
nantly defined by the components close to the k-space 
origin. Objects or organs with sharp high-contrast edges 
contain a significant spectral density at the k-space periph- 
ery. The majority of biological samples show very local spec- 
tral density in k-space, centered around k= 0. Some organs 
and tissues, such as brain cortex, have a fractal-like nature“ 
and therefore show a slower decay of spectral density in 
k-space. 

To understand the properties of k-space and its rela- 
tion to image/object space, it is important to develop an 
intuition regarding the local and global properties of these 
two spaces. This topic is reviewed in detail by Paschal and 
Morris.” K-space describes the object using a set of global 
planar waves. Every sample in k-space describes the contri- 
bution of the wave with a corresponding frequency to the 
entire image. Therefore, a change in a single sample in 
k-space affects theoretically the entire image. Similarly, a 
change in the intensity of a single pixel (e.g., in a dynamic 
process) generally affects all k-space samples. However, when 
such change can be allocated to a larger cluster of pixels, 
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predominantly the center of k-space is affected. The latter 
phenomenon is the basis of a range of data-sharing techni- 
ques, which achieve a higher apparent temporal resolution 
in dynamic imaging by updating the central region of 
k-space more frequently. ? 

Different strategies to populate k-space with measured 
data are termed “k-space sampling trajectories” or “k-space 
trajectories”. The most common and clinically relevant 
approach collects data on a rectilinear grid in k-space (so-called 
“Cartesian” sampling), because it allows for a computationally 
efficient image reconstruction using the fast Fourier transform 
(FFT). Even though we primarily consider the Cartesian sam- 
pling strategy, it is necessary to distinguish between the multi- 
tude of k-space trajectories, which differ from each other in the 
way and the order the nodes of the rectilinear 2D or 3D 
k-space grids are visited. Other popular data acquisition 
schemes include radial,’ PROPELLER, "! and spirals.'* 

Typical motion-induced deterioration effects observed 
in MR images consist of a combination of the following 
basic effects: blurring of sharp contrast or object edges 
(intuitively similar to photography), ghosting (both coherent 
and incoherent) originating from moving structures, signal 
loss due to spin dephasing or undesired magnetization evo- 
lution, and the appearance of undesired strong signals. 

The first two points are related to the signal readout 
process, whereas the latter two are related to the signal gen- 
eration and contrast preparation within the pulse sequence. 
Here we refer to ghosting as a partial or complete replica- 
tion of the object or structure along the phase-encoding 
dimension, or along multiple phase-encoding dimensions 
for 3D imaging. Periodic motion synchronized with the 
k-space acquisition results in a coherent ghosting with the 
number of replicas corresponding to the frequency of & 
space modulation: two ghosts result if every second line is 
altered, four if every fourth, etc. Deviations from perfect 
periodicity in &-space result in incoherent ghosting, appear- 
ing as multiple overlapped replicas, and sometimes seen as 
stripes in the phase-encoding dimension. 

The main cause of readout-related motion artifacts is 
the inconsistency between the various portions of the & 
space data used for the image reconstruction or between the 
data and the signal model assumed in the reconstruction. 
Simple reconstruction using an inverse FFT (iFFT) assumes 
the object has remained stationary during the time the & 
space data were sampled. A violation of this assumption 
results in artifacts. 

To demonstrate possible motion artifact appearance we 
have simulated a range of motions and k-space reordering 
schemes. Figure la shows the digital phantom used based 
on the Shepp-Logan phantom with altered intensities and a 
grid of small objects added. Figures 1b-d demonstrate the 
range of simulated motions (continuous rotation to 10°, 


continuous translations by 10 pixels in vertical and 
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FIGURE 1: (a) Original FFT reconstruction without motion, (b-d) sum of first and last source images to show the range of simu- 
lated motions for (b) continuous monotonic rotation to the total angle of 10°, (c) continuous vertical translation, and (d) continu- 
ous horizontal translation (in both cases 10 pixels for a matrix of 256 x 256). Periodic horizontal translation has the same 
amplitude and is not shown. Images (e-h) demonstrate the results of averaging the motion as would correspond to photography 
with an equivalently long exposure time; (e) rotation, (f) vertical translation, (g) horizontal translation, and (h) periodic horizontal 
translation. Note the loss of detail and edge information due to the motion blur and the enhanced blurring in (h) due to the fact 
that with sinusoidal motion the object spends more time close to the terminal positions. Images (i-l) show simulated MRI acquisi- 
tions with linear k-space ordering. Images (m-p) show simulated two-shot interleaved MRI acquisitions. Plots of the motion as a 


function of k-space position for the representative images are presented in Figure 2. 


horizontal directions (Figs. 1b-d, respectively) and a peri- 
odic sinusoidal oscillation in the horizontal direction with 
the same amplitude (not shown). To provide a reference of 
image quality, a photography-like image acquisition has 
been simulated, averaging the motion in image space (Figs. 
le-h). It is noteworthy that the image lh appears more 
blurred than 1g, because for periodic motion the object 
spends more time in proximity of the terminal positions. 

It is common knowledge in the imaging community 
that MRI is particularly sensitive to motion. However, in 
the case of slow continuous drifts and sequential k-space 
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ordering, Fourier acquisition is actually significantly less 
motion-sensitive than a corresponding “photographic” image 
taken with the same “exposure time” (Figs. li-k). Periodic 
motion however produces very strong ghosting (Fig. 11). In 
the case of interleaved multishot k-space ordering, even slow 
continuous drifts produce significant ghosting (Figs. 1m—o). 
Note that ghosting reduces for the periodic motion in 
Figure lp in comparison to Figure 11 because the inter- 
leaved acquisition scheme traverses k-space faster, which 
reduces the number of full oscillations effectively seen in the 
complete k-space data set. If the interleaving scheme is 
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effective x-translation 


translation 


k-space position 


effective x-translation 


k-space position 


effective x-translation 


translation 


b k-space position 


effective x-translation 


translation 


d k-space position 


FIGURE 2: Effective translation (as a fraction of the FOV) as a function of the k-space position for (a) Figure 1j and 1k, (b) Figure 
1l, (c) Figure 1n and 10, and (d) Figure 1p. As seen, if the external disturbance oscillates as a function of k-space position, e g. in 
cases (b), (c), and (d), ghosting will result (see corresponding images in Fig. 1). Reduced ghosting in Figure 1p can be explained 
by the effective reduced frequency of oscillation of the position as a function of the k-space coordinate (compare plots b and d). 


Plots of rotations are omitted due to their similar appearance. 


synchronized with the underlying periodic motion, the arti- 
facts reduce further to the level comparable with that in Fig- 
ure 1k. This phenomenon is often exploited in cardiac 
imaging and is discussed in more detail below. Plots of the 
motion as a function of k-space position for the representa- 
tive images in Figure 1 are presented in Figure 2. 

Although slow continuous motions do occur in clin- 
ical settings, for example due to the gradual relaxation of 
the neck muscles in head imaging, they only strongly 
affect sequences relying on interleaved k-space acquisitions, 
e.g., T2-weigted TSE/FSE (Turbo/Fast Spin Echo) imag- 
ing. In addition to periodic processes (e.g., breathing, car- 
diac motion, blood pulsation and tremors), sudden 
position changes, for example due to swallowing in head 
imaging or insufficient breath hold capability in abdomi- 
nal scanning are the effects that often lead to artifacts. 
The effect of sudden motion on the Fourier acquisition is 
visualized in Figure 3. Here, single 10° rotation events 
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have been simulated, with the motion occurring at differ- 
ent time points, leading to varying proportions of k-space 
data that are inconsistent. Figures 3a—d and Figures 3e—h 
assume linear and interleaved k-space ordering, respec- 
tively. In both cases 12.5% (1/8th) of the inconsistent k- 
space data result in a negligible effect; however, image 
quality deteriorates much faster for the interleaved acquisi- 
tion. Centric k-space reordering demonstrates even higher 
resistance against artifacts (Figs. 3i-I). 

The following conclusions can be drawn from these 
simulations: (i) MRI with a sequential k-space ordering is 
fairly insensitive to slow continuous motion; (ii) motion as 
a function of k-space position is the most relevant parameter 
for the appearance of artifacts; (iii) it is important where in 
k-space fast motion occurs, as data corruption near the cen- 
ter of k-space produces stronger artifacts than data corrup- 
tion near the k-space periphery; (iv) repeated modulation in 
k-space amplifies artifacts and should be avoided. 
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FIGURE 3: Simulations of a single sudden orientation change during the k-space acquisition for different k-space acquisition strat- 
egies and different amounts of inconsistent k-space data. For a linear k-space ordering, images reconstructed from datasets con- 
taining 12.5% (a), 25% (b), 37.5% (c), and 50% inconsistent data (d) have been simulated. Hardly any artifacts are visible if the 
corruption occurs less than half of the maximum k-space distance to the center. Similar conclusions follow for the two-shot inter- 
leaved k-space acquisitions with 12.5% (e), 25% (f), 37.5% (g), and 50% of inconsistent data (h). i-l: Images present a similar simu- 
lation for the centric reordering with (i) 50% of data inconsistent with respect to the k-space center data, (j) 62.5%, (k) 75%, and 
(I) 87.5%, respectively. The precise threshold at which visually detectable artifacts start to appear strongly depends on the image, 
primarily on the presence of small high-contrast features, but the overall tendency is that below a certain distance from the k- 


space center a single motion event produces only negligible artifacts. 


Types of Motion 

The different types of motion relevant for in vivo MRI can 
be loosely separated into the following categories: rigid body 
motion, elastic motion and flow. 

Rigid body motion, or “bulk motion”, includes 1D 
translations (requiring a single parameter for its mathematical 
representation), multi-dimensional translations (requiring sev- 
eral parameters), and completely unconstrained rigid motion 
(requiring six parameters). For example, diaphragm motion is 
sometimes approximated using a one-dimensional translation, 
but involuntary head motion normally requires six degrees of 
freedom, comprising three translations and three rotations. 

Elastic motion typically includes stretching, compres- 
sion and shearing along three axes, in addition to rigid body 
motion. It requires 12 degrees of freedom for its complete 
representation. Elastic motion is observed in the abdomen, 
where various locations experience different displacements 


and deformations. Flow can in some cases be assumed to be 
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one-dimensional, for example for laminar blood flow in 
small or medium-sized vessels or CSF flow in the cervical 
spine. It then requires only a single parameter (velocity) 
for its complete representation. More complex flow is 
observed in the heart and aorta. This requires a two- or 
three-dimensional velocity vector field for its complete 
description. 

The most complete description is provided by the dis- 
placement field, which consists of the trajectories of each 
particle within the object as a function of time. Although 
being complete, the displacement field is often not measura- 
ble with sufficient precision. Luckily, for the vast majority 
of practical situations this description is too detailed. Physi- 
cal, mechanical and physiologic constraints placed on the 
human body (e.g., incompressibility of liquids and a major- 
ity of tissues, presence of bones and joints, etc.) allow for a 
drastic reduction of the number of parameters required for 
the description of motion. 
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FIGURE 4: Velocity effects on the phase of MR signals. Stationary spins can be refocused with a pair of gradient pulses of equal 
area and opposite polarity, which is a basis of gradient echo. Spins moving along the gradient direction acquire additional phase 
due to the fact that dephasing and refocusing occur at different physical locations where the same gradient indices different fre- 
quency shifts, leading to different amounts of the phase acquired during dephasing and refocusing periods. Greyed segments on 
the corresponding phase circles mark the amount of phase acquired by the moving spins during the last gradient pulse. 


Data Corruption Mechanisms 
Several physical mechanisms contribute to the appearance of 


motion artifacts. The most important and relevant for in 


vivo MRI are the following: 


i 
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Incorrect phase accumulation due to tissue motion dur- 
ing the periods in the pulse sequence where gradients 
are switched on. MRI relies on the ability to create gra- 
dient echoes, and typically assumes that tissues are sta- 


ii. 


tionary during imaging. As seen in Figure 4, spins 
moving in the direction of the gradient acquire addi- 
tional phase. If this phase varies for different phase- 
encoding steps, inconsistencies in k-space arise; 

The so-called excitation history effect appears when 
slice-selective RF pulses are used for excitation, satura- 
tion or refocusing. As seen in Figure 5, out-of-plane 
affects the desired 


motion between such pulses 
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time 


FIGURE 5: Excitation history effects visualized for a three-pulse MR pulse sequence. In the case of no motion, regions excited by 
all three pulses overlap resulting in the desired signal evolution. The respective brain regions are marked using half-transparent 
bands with the color of the corresponding slice-selective RF pulse. Overlap of red and green results in light brown. When light 
brown is followed by an overlap with blue, a grey color results, which indicates the desired spin evolution. In the case of motion, 
(rotation after the first pulse, followed by a combination of rotation and translation after the second pulse, original position is 
shown as dotted line) the grey region corresponding to the desired spin evolution is substantially reduced, while colored regions 
with undesired signal evolution appear. Therefore, desired signals are decreased, and undesired signals are produced, resulting in 
artifacts, such as ghosting or signal disturbances. The unwanted signals may also cancel some of the desired signals, causing signal 


voids. 


evolution of the signal and typically results in genera- 
tion of signals that are too strong or too weak, or causes 
alternations in image contrast. For multi-shot imaging 
this results in magnitude inconsistencies in k-space; 

iii. Motion of the tissues and body parts often affects mag- 
netic fields used in MRI. The BO field within the imag- 


ing volume changes significantly if the orientation of 
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a b 


FIGURE 6: The effect of rotations on the k-space sample loca- 
tions in object coordinates. In the absence of rotations (a), the 
k-space samples are arranged on a uniform grid fulfilling the 
Nyquist criterion. Rotation of the object during the image 
acquisition rotates the spatial frequency components associ- 
ated with the object with respect to the encoding gradients, 
which effectively redistributes the samples in k-space (b). This 
leads the sampling density in some areas of k-space to fall 
below the Nyquist limit and results in streaking and ghosting 
artifacts in the images, which are difficult to correct for in 
reconstruction even if the motion is known. 
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the body part changes but may also be affected by 
motion of relatively distant body parts due to long- 
ranging magnetic susceptibility effects.'*'4 Both trans- 
mit and receive B1 fields may also change depending 
on the body position or due to the motion of the sur- 
face coils themselves. 

iv. Rotations during the k-space data acquisition process in 
multi-shot imaging produces inconsistencies between 
the effective readout directions for different k-space 
lines (see Fig. 6). Even if the true motion during the 
acquisition is known and compensated for, the resulting 
k-space is not sampled homogeneously, which may 
result in ghosting or streaking artifacts. 

v. Cardiac or respiratory motion leads to so-called 
“physiologic noise” in fMRI. Because the fMRI signal is 
very small," any other dynamic changes caused by, for 
example, periodic blood flow, appears as a confounding 
factor in the analysis. The situation is further exacer- 
bated by the relatively low temporal resolution of typi- 
cal fMRI acquisitions leading to a spectral aliasing of 
the physiologic signals. 


Data inconsistencies between k-space lines or segments 
in multi-shot imaging typically result in ghosting or blur- 
ring. Most of the above mechanisms may also affect the 
contrast preparation modules (e.g., diffusion weighting or 
inversion) and thus generate undesired contrast variations 
without visible ghosting artifacts. This occurs particularly in 
single-shot imaging. 
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TABLE 1. Common notion artifact mitigation strategies 


Motion prevention 
Training 

Distraction 

Feed and wrap (for babies) 


Foam restraints 


Artifact reduction 

Faster imaging 

Insensitive sequences 
Gradient moment nulling 


Saturation bands 


Sedation Triggering and gating 
Bitebars/head holders Phase reordering 
Breathhold 


Motion correction 
Navigators 
Self-navigated trajectories 
Prospective correction 


Retrospective correction 


Artifact Mitigation Strategies 


As discussed above, there are many types of motion and 
many possible mechanisms for data corruption resulting 
from each motion type. It follows, therefore, that many dif- 
ferent strategies have been developed to address these chal- 
lenges. In the following, we categorize these methods into 
three distinct groups: motion prevention, artifact reduction 
and motion correction (Table 1). 

Motion prevention is the most obvious method of 
suppressing motion artifacts. If motion can be avoided, then 
the effects discussed early in the review are prevented and 
other, more complex, strategies become redundant. Unfortu- 
nately, preventing motion is not always practical and so arti- 
fact reduction or motion correction strategies must often be 
used. Artifact reduction is defined in this review as acquisi- 
tion strategies that reduce artifacts in the resulting image or 
replace them with those with a less dramatic appearance, 
when compared with standard Cartesian acquisition 
methods. Motion correction, on the other hand, normally 
involves the explicit estimation and compensation of 
motion. Combinations of all three methods are possible, 
just as multiple tools may be used together in a complemen- 
tary manner in the toolbox analogy mentioned above. 


Motion Prevention 

Prevention is the most commonly attempted approach to 
avoid motion artifacts in clinical MRI. The main prevention 
tools available are briefly summarized here. 

For young children, training with a mock MRI is useful 
to avoid bulk motion by reducing anxiety. However, this 
method is not widely available and is not particularly cost 
effective. In infants, imaging can be timed after feeding to 
take advantage of sleeping. Babies can be well wrapped, which 
also encourages sleeping. Foam restraints or special inflatable 
devices are in common use, the latter particularly in the case 
of infants. However, these methods are only partially effective. 
Nevertheless, Windram et al. have recently reported failure- 
free imaging of 20 young patients with a congenital heart dis- 


ease using the feed-and-sleep technique.'® 
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Moderate sedation, or in some cases, general anesthesia, 
is often used '”’'® when imaging young children. Surprisingly 
enough, sedation is not always effective, with occurrence of 
motion artifacts in approximately 12% of cases '’ and an 
increased prevalence of failures in older children. According 
to Malviya et al.,'” general anesthesia is more reliable, as it 
has a much lower incidence of motion artifacts (approxi- 
mately 0.7%) but is associated with much greater costs (in 
part due to the requirement of a presence of an anesthetist), 
more severe side effects and greater health risks. 

In the early days of MRI, bite bars mounted on head 


© . . . 
19.20 were deemed an effective means of avoiding head 


coils 
motion. However, due to their cumbersome set up and sig- 
nificant discomfort, these approaches have not found wide 
acceptance in clinical practice. The same applies to other 
immobilizing systems, such as a personalized plaster cast 
head holder.2! Vacuum immobilizer cushions, which are 
available from several suppliers, address the set-up issues, 
but, in the experience of the authors, are perceived as very 
uncomfortable after a time period of 10 to 30 min. 

In the case of abdomen or cardiac imaging, breath- 
holding is an obvious, and effective, methods of suppressing 
motion artifacts due to breathing (see Chapter 12.3 in Bern- 
stein et al.).°? Breath-holding, however, restricts scan dura- 
tion severely (between 10 and 20 s in clinical settings), 
which in turn limits the image quality, resolution and. cover- 
age. Image quality is further restricted by the quality of 
breath-hold (e.g., presence of drifts or tremor). To overcome 
time limitations multiple breath-holds are often used despite 
limitations associated with breathing hysteresis and limited 


reproducibility of breath-hold positions. 


Artifact Reduction 

The basic remedial approaches and tools to reduce motion 
artifacts were devised approximately 30 years ago, in the 
early days of MRI. A 1986 review article by Bellon et al. 7° 
lists strategies such as shortening imaging time, optimizing 
the phase encoding direction and ordering, and reducing 
flow sensitivity. These basic approaches remain valid to date, 


Volume 42, No. 4 


ASUDOT] SUOWUOJ ANLAN aquorfdde oy Aq pouJdAos are sAN WO ‘asn Jo SIMI 107 AreAqr] IUUQ ÁLMA UO (SUONIPUOD-puL-sUIOY/WIOD ÁM AIeIgIaUTUo//:sdyy) SUONIPUOD pUe SUIL A 99g ‘[EZOZ/TI/LTZ] VO ArerQr] IUUQ Ad|IAA ‘VIsoUOpUT AOIg WN Aq OSSPT HUL/ZOOT O [/HOp/wios Aoi Arerqrouruoy/:sduy wos popeopumog ‘fp ‘STOT ‘YSSZZZST 


but new ways of achieving the same goals have appeared. 
Therefore, for the sake of completeness, we review the cur- 
rent artifact reduction strategies below. 


FASTER IMAGING. In analogy to photography, shortening 
the image acquisition time below the typical time scale of 
motion is the most straightforward approach to reduce 
motion artifacts. The wider availability of fast imaging techni- 
ques is largely responsible for the improvements seen in MRI 
with regard to the immunity to motion, rather than the more 
complicated motion correction techniques discussed below. In 
the case of spontaneous movements, shorter scan times mean 
that the patient is less likely to become uncomfortable and 
move. In the case of cardiac or breathing motion, fast sequen- 
ces acquire more data within a single period. An early example 
of fast gradient echo imaging is the FLASH (fast low-angle 
shot) sequence,” which uses low-angle excitation and spoils 
transverse magnetization between TRs. Progress in gradient 
hardware has resulted in the wide availability of fully- 
refocused gradient echo sequences 25 (so-called balanced 
steady-state free precession, bSSFP). For the same reason, 
echo planar imaging (EPI), which was initially deemed 
barely feasible, has became a workhorse in functional neuroi- 
maging. T2-weighed imaging has became clinically feasible 
due to the adoption of the echo-train principles,” but still 
remains too slow for certain applications, such as breathhold 
abdominal imaging. In such cases a single-shot fast spin-echo 
(a.k.a. HASTE) is often used as a means of freezing motion at 
expense of a moderate loss in spatial resolution. 

Parallel imaging’ ° has revolutionized many MRI 
applications by shortening acquisition times, typically by a 
factor of 2 to 3. More recently, these approaches have been 
extended to acceleration in two encoding dimensions??? lead- 
ing to even higher gains for 3D imaging or novel approaches, 
such as simultaneous multi-slice imaging.” The introduction 
of compressed sensing MRI * in 2007 has provided a further 
tool for acceleration of MRI acquisitions and has paved a way 
to new motion-compensated applications.” 

Fast imaging, however, compared with the standard 
counterpart, frequently represents a compromise between 
the acquisition time and other confounding issues. Single- 
shot imaging typically has a limited spatial resolution. EPI 
suffers from geometric distortions and chemical-shift arti- 
facts. Parallel imaging methods rely on the accuracy of the 
coil calibration data and may result in new types of artifacts 
if motion happens during the acquisition of the reference 
data.*°?* Compressed sensing is often criticized for the 
introduction of the new types of artifacts (e.g., patchiness, 
synthetic appearance), which may mask important anatomic 
features in some cases.” 

GRADIENT MOMENT NULLING. Gradient moment nul- 
ling is a widely used method to reduce artifacts produced by 
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FIGURE 7: The effect of linear motion on the phase acquired 
by the spins under the action of the encoding gradients (fre- 
quency encoding in this example). A bipolar gradient (a) refo- 
cuses the stationary spins at the moment when the net area 
under the gradient reaches zero. This is because the phase 
accrual in this case is proportional to a product of the gradient 
amplitude and the time duration. Contrary to that, the phase 
acquired by a spin moving along the gradient direction is pro- 
portional to the gradient amplitude, duration and the displace- 
ment of the spin. Because the displacement under continuous 
motion is proportional to the time, the phase shows a quad- 
ratic dependency on time. Therefore, a bipolar gradient fails to 
refocus moving spins. Consequently, such spins will be poorly 
located and a motion artifact will be seen in vascular structures 
and moving fluids. It is possible to account for the quadratic 
behavior of the phase by introducing a third gradient lobe, as 
shown in (b). Here, to keep the polarity of the frequency 
encoding gradient the sign of the two preceding lobes has 
been adjusted. The gradient area ratio of 1:-2:1 does not 
induce additional dephasing for the stationary spins and allows 
for refocusing the spins moving with a constant velocity. Other 
gradient schemes are possible taking into account gradient 
amplitude and slew rate limits, but velocity compensation 
always increases the minimal echo time. 


objects or tissues moving with a near-constant velocity. It is 
particularly useful for the suppression of flow-related artifacts. 
The method attempts to minimize the phase shifts acquired 
by both stationary and moving spins at the echo time. With- 
out gradient moment nulling (Fig. 7a), stationary spins are 
refocused, but spins moving along the direction of the applied 
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gradient are not. With gradient moment nulling (Fig. 7b), 
both stationary spins, and spins moving with a constant veloc- 
ity, are refocused. Mathematically this is equivalent to mini- 
mizing both the Oth gradient moment (area under the 
gradient waveform) and the Ist gradient moment at the echo 
time. The gradient moment nulling condition is exactly or 
approximately fulfilled for a range of multi-echo imaging 
sequences *° such as multi-gradient echo or turbo spin echo, 
but only for even echoes within the echo train. This phenom- 
enon is termed “even echo refocusing” in the literature. 
Incorporation of gradient moment nulling almost 
always slows down the sequence as it requires additional gra- 
dient lobes to be added. Currently the majority of scanner 
vendors provide this feature for one or several encoding 
dimensions as an option amongst the sequence parameters, 


typically referred to as “flow compensation”. 


MOTION INSENSITIVE SEQUENCES. Imaging methods 
relying on spiral or radial encoding strategies are often more 
motion-tolerant than Cartesian sequences, in that the 
motion artifacts have a more benign appearance in the 
images.*' The reason for this is that these strategies typically 
oversample the center of k-space, which either smears out 


or may be 
5-49 


the artifacts, due to the inherent averaging “1~“4 
used for more advanced correction in postprocessing.“ 
Furthermore, if the readout starts at or close to the center 
of k-space, it is easier to achieve gradient moment nulling, 
which reduces phase errors between the shots. 


PHASE REORDERING. The order of acquisition of k-space 
lines can make a substantial difference to the appearance and 
severity of motion artifacts when imaging a moving object. 
Changing the ordering to minimize artifacts is particularly 
common to prevent artifacts from respiratory motion. An 
early example of this is respiratory-ordered phase encoding 
(ROPE), where neighboring lines in k-space are acquired at 
points in the respiratory waveform that are also close 
together. This minimizes side lobe peaks in the point-spread 
function (PSF), which reduces ghosting. Such techniques are 
often generally referred to as “respiratory compensation”. A 
major advantage of respiratory compensation over respiratory 
gating is a much shorter scan time. 

Respiratory compensation requires information about 
motion caused by the patient breathing. Typically, respira- 
tory bellows are used to obtain a waveform that gives infor- 
mation about the patient’s chest position by means of 
measurement of the movement of air in the bellows. It is 
important that the bellows are correctly attached to measure 
the full range of respiratory motion. MR navigators that 
track breathing motion are also often used as a source of 
information for the phase reordering.”! 


SPATIAL SATURATION BANDS. Movement of parts of 
the anatomy (e.g., blood flow from arteries, breathing 
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motion in abdominal imaging, etc.) can generate motion 
artifacts that overlay the organ of interest.°”°? Spatial satu- 
ration bands can be used to suppress signal from such mov- 
ing tissue. To achieve this, one or more spatially selective 
90° RF pulses are applied before the actual excitation pulse. 
The transverse magnetization is then dephased by spoiler 
gradients. Moving spins from the saturated region are, there- 
fore, prevented from contributing to artifacts in the imaged 
volume. However, the additional RF pulses lead to an 
increase in SAR and measurement time — clinically, this 
results in a decreased number of slices for a given TR. 


TRIGGERING AND GATING. Patient motion resulting 
from physiologic effects such as heartbeat or breathing is 
generally periodic. This is an advantage for motion correc- 
tion, as MR data acquisition can be timed to the different 
points within the cardiac cycle ** or breathing position.” 

There are two general approaches to enforce consis- 
tency of the data acquired in presence of (quasi-) periodic 
motion. Data can either be collected at the same point in 
the cycle (triggering) or acquired continuously and then 
reordered retrospectively (gating). Triggering is often easier 
to implement and it has the advantage of a more precise 
synchronization with the underlying motion. On the other 
hand triggering by necessity introduces a delay, which both 
disrupts the signal steady state and misses some part of the 
dynamic cycle in applications like cardiac cine imaging. 
Gating is free of these shortcomings, but it has to sacrifice 
some scanning efficiency and allow for a certain degree of 
redundancy to ensure that sufficient data are available for 
retrospective image reconstruction. 

Both triggering and gating rely on an additional signal 
or several signals correlated to the physiologic motion in 
question. For cardiac imaging, electrodes are normally 
placed on the subject’s chest to record the electrocardiogram 
(ECG). This has the drawback of slower patient workflow 
and has some safety implications.*° At higher field strengths 
(3 T and above) the quality of the ECG trace may be 
reduced substantially, so it is, therefore, important to opti- 
mize the placement of the EEC electrodes and leads.”’ 
Another possibility is to use a pulse oximeter, typically 
attached to subject’s finger. However, peripheral signals are 
less sharp and precise and are more suitable for gating 
because of the delay associated with a limited velocity of the 
pulse wave in the human body. Respiratory gating and trig- 
gering can rely on respiration sensors, such as belts or bel- 
lows, typically affixed to the subject’s chest or belly, as 
discussed above. However, the reliability of such devices and 
the relation to the breathing-induced motion of the internal 
abdomen organs are substantially lower than in the case of 
cardiac motion. Therefore, MR navigators following the 
position of the diaphragm or the organ of interest are often 
preferred. Navigators are discussed in more detail in a 


Volume 42, No. 4 


ASUDOT] SUOWUOJ ANLAN aquorfdde əy) Aq pouJdAos ase sAN WO ‘asn Jo SIMI 107 AreQr] IUUQ Adi AA UO (SUOTTPUOD-puL-SULIOY/WIOD ÁM AIeIgaUTUO//:sdyy) SUONPUOD pUe SUIL A 99g ‘[EZOZ/TI/LT] UO ArerQry] IUUQ Ad|IAA *LISIUOpPUJ AOIg IN Aq OSSPT UWUL/ZOOT O 1/10Pp/W02 Kom ArerqrouTuoy/:sduy Wor popeopumogd “p ‘STOT ‘YSSZTZZST 


separate section below. Another popular approach acquires 
one or several nonencoded samples following the RF excita- 
tion. This is commonly referred to as FID-gating or self- 
gating. >? 

Generic limitations to both triggering and gating arise 
from the violation of the two basic assumptions behind 
these techniques: periodicity of the disturbing motion and 
stationarity of the imaged object or organ. Therefore, car- 
diac MR imaging in patients with arrhythmia remains chal- 
lenging. In pediatrics both breathing and heartbeat may be 
very irregular, which limits the achievable image quality 
severely. Also in adults free-breathing exams fail sporadically, 
due to breathing position drifts and changes in breathing 
patterns. Imaging of fast nonstationary and irreproducible 
events, such those observed in cardiac perfusion, is an area 
of active research (e.g., Schmidt et aij 


Motion Correction 


MR NAVIGATORS. MR navigator methods acquire the 
required position information during an imaging sequence 
using the MR scanner itself.°'“®? This is achieved by fre- 
quently playing out ultra-short (on the order of millisec- 
onds) navigator sequences during the imaging process. 
These navigator sequences allow for determination of the 
object’s position by comparing each acquired navigator sig- 
nal to a reference. 

Generally, navigators can be separated into two catego- 
ries: those that work in k-space and those that work in 
image-space. They can be further classified in terms of their 
dimensionality: 1D, 2D, or 3D. Because any additional 
information encoded by the navigator comes at a cost of 
increased acquisition and processing times, and increased 
complexity, navigators are typically simplified to an extent 
acceptable for the target application. Therefore, 1D and 2D 
navigators are popular especially in abdominal and cardiac 
imaging to correct for breathing motion, by tracking the 
position of the diaphragm.°* In brain imaging, 2D navi- 
gators allow for a correction of rigid body motion in the 
navigator plane.°* ©” In DWI even the smallest movements 
caused by brain pulsation can lead to phase changes due to 
the strong gradients used for diffusion weighting. Variations 
of the image phase leads to substantial image artifacts in 
multishot-DWI. MR-navigators measuring the image phase 
after diffusion weighting (in 2D) allow the correction of 
these phase instabilities.”°7' 

3D navigators allow for the correction of patient 
motion in up to 6 degrees of freedom.°*”*”* They also 
require substantial time to be available in the host sequence, 
on the order of hundreds of milliseconds. This time is avail- 
able in many sequences, due to contrast preparation or 
recovery delays; therefore, 3D navigators have become 
increasingly popular, especially in brain structural imaging. 
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Spectroscopic acquisitions offer even more “empty space” 
within the sequence, which is sufficient to acquire several 
navigator echoes, or even specialized shim navigators to cor- 
rect for shim alterations due to motion.’””° Recently, sub- 
stantially faster 3D navigators, relying on parallel imaging 
or multiband excitation,” have been presented at various 
conferences and workshops. 

To improve the accuracy of the navigation sequence, 
and to reduce interference with the imaging process, it is 
possible to use additional spatial-frequency-tuned markers,” 


POT oes signals of the 


miniature radio-frequency probes, 
endogenous fat.°”*' However, many of these methods were 
so far restricted to 1D ®© and 2D,’”*! and are generally 
only used in research applications. 

An alternative to the navigator methodology is self- 
navigated methods, which allow the determination of motion 
from the MR imaging data themselves.“ 74>58_A]] of these 
methods are very specific to a particular imaging scheme and 
are often limited to 2D acquisitions. One popular variant of 
self-navigated sequences is PROPELLER, which is discussed 
in detail in the next section. Self-navigation methods can also 
be extended to 3D, for example in conjunction with a three 
dimensional radial imaging scheme.*® 

Navigator information can be used either retrospec- 
tively to correct the acquired imaging data 8° or prospec- 
tively to adapt the imaging scan to the patient’s position. 
Note that for prospective motion correction, motion param- 
eters need to be extracted from the navigator data in real 
time. Advantages and disadvantages of both prospective and 
retrospective approaches are discussed below. 

A detailed overview of different navigator approaches 
can be found in Chapter 12.2 of Bernstein et al.’” 


PROPELLER 


The acronym “PROPELLER” is an abbreviation for 
“periodically rotated overlapping parallel lines with 
enhanced reconstruction”.'' PROPELLER imaging, and its 
variants, follow a strategy similar to that of radial imaging. 
However, instead of acquiring the data projection-by- 
projection, strips of several parallel k-space lines are col- 
lected. This can be realized using different readout strategies 
like RARE,*? EPI,°*®° or GRASE.™* This strip is rotated 
around the k-space origin in the subsequent acquisition 
cycles until a disc of k-space is filled. 

This sampling scheme provides substantially reduced 
motion sensitivity due to the strong oversampling of the & 
space center and its radial character. In addition, the robust- 
ness to motion can be further improved during the recon- 
struction.*° Low-resolution images reconstructed from single 
blades can be used to quantify in-plane translations and cor- 
rect for the resulting inconsistencies between the blades before 
the final image reconstruction. Additionally, the central disc 
of &-space—which is acquired by all strips—can be used to 
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detect and correct in plane rotations of the single k-space 
strips. By cross-correlating the central disc of the corrected 
strips, through-plane motion can be detected. Such motion 
cannot be corrected for retrospectively; however, by introduc- 
ing a factor that puts less weight on corrupted blades, the 
influence of through-plane motion can be reduced. 

The PROPELLER sampling strategy can be combined 
with undersampling techniques.*” It can also be used to intro- 
duce diffusion weighting to the RARE readout module,*® 
which was shown to be of interest for DTI at high fields.®° 

The PROPELLER methodology has been shown to be 
a very useful tool in clinical settings °° and is available on 
the majority of imaging platforms (under different names). 
Drawbacks of PROPELLER are the increased image acquisi- 
tion time, due to the strong oversampling around the & 
space center, the limitation to 2D imaging and a limited 
robustness against through-plane motions. 


Prospective Correction 

Prospective motion correction is an intuitive approach where 
the relative spatial position and orientation between scanning 
coordinates and the object of interest is kept constant. For the 
correction of head motion, a rigid object is assumed, moving 
in 6 degrees-of-freedom. If the position and orientation of 
such an object can be measured in real time, the magnetic gra- 
dients, RF pulses, and receiver frequency and phase can be 
adjusted accordingly. Rotations of the object require a rotation 
of the encoding gradients; translations require a change of 
transmit and receive frequency and phase.”! 

In the case of head motion, head pose can be meas- 
ured using MR navigators, working either in image space 
(e.g., PROMO 92:93 PACE °$), or k-space (e.g., cloverleaf, 
spherical or orbital navigators°°°”’”), Alternatively, an exter- 
nal tracking device can be used, including stereo camera sys- 
tems,” miniature RF probes, ’*”??° in-bore camera 


systems, RR 


or ultrasound systems.'°? Navigators need to 
be compatible with the sequence timing and typically have a 
low temporal sample rate. On the other hand, external 
tracking typically requires an MR compatible design and 
accurate calibration to the scanner coordinate system. Fur- 
thermore, the majority of the external approaches require a 
tracking device to be mounted on the subject, which is a 
concern for routine applications. 

The topic of prospective motion correction has gained 
popularity in the last few years resulting in numerous new 
applications including fMRI, '°™>!°? DWI, "0104 and spec- 
troscopy. °% !°7 Although it is an extremely promising 
approach for neuroimaging, it does have some limitations, 
including practical considerations (e.g., marker fixation for 
external tracking systems) and uncorrectable effects (e.g., 
motion-related BO distortions).'°% Covering this topic in 
detail is beyond the scope of this study; interested readers 


i . 4109 
are referred to a recent review article. 
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Postprocessing Techniques and Retrospective 
Correction 

While prospective correction methods aim at maintaining 
data quality during the acquisition, retrospective techniques 
attempt to improve data consistency afterward, by modify- 
ing the collected data or the reconstruction model. This can 
be done by either including the exact knowledge about the 
motion during the scan (e.g., using navigator data) or with 
iterative algorithms (e.g., optimizing image entropy, gradient 
entropy or other measures of artifacts).''0-'! 

The basic idea behind these methods is to undo the 
motion-related changes that occurred to the MR data. For 
rigid body motion these are described according to the Fou- 
rier theorems: a translation of the object leads to a phase ramp 
in the acquired k-space, an object rotation corresponds to a 
rotation of &-space.''* While translations are relatively easy to 
correct by applying a phase change to the acquired data, the 
correction of rotations requires the use of non-Cartesian 


; 114,11 
reconstruction methods ° 


116-119 


and includes some sophisti- 
cated algorithms, which are computationally intensive. 

In 3D imaging, arbitrary rigid body motion can be cor- 
rected, as long as the acquired signals are only damaged but 
not lost (e.g., there are no signal dropouts due to intra-voxel 
dephasing). For 2D imaging these correction methods are 
limited to in-plane motion, as through-plane motion during 
the scan results in inconsistencies, which cannot be corrected. 

Elastic motion remains a major challenge for both pro- 
spective and retrospective correction approaches. Especially 
in abdominal and cardiac MRI, the complexity of the 
underlying motion restricts the current concepts to a combi- 
nation of gating/triggering followed by a correction of affine 
motion within the gating window.°”'”° 

An alternative approach to retrospectively improve the 
quality of datasets which are only partially affected by 


1.1?! By discarding 


motion is described by Bydder et a 
motion-corrupted data and filling the resulting parts of & 
space using parallel imaging techniques, the quality of the 
reconstructed images could be increased at a cost of SNR 
and possibly some residual blurring. 

By including a general description of motion into the 
MR signal equation’ using measured or estimated dis- 
placement fields’? free-breathing cardiac imaging has been 
demonstrated. Recently this approach has been extended to 
incorporate compressed sensing methodology.” Unfortu- 
nately such retrospective image reconstruction methods are 
extremely computationally intensive, which presently 


restricts their clinical acceptance. 


Conclusions 

Due to peculiarities of the image acquisition process in 
MRI, motion causes a varying range of artifacts, including 
blurring, ghosting, signal dropouts and unwanted signal 
enhancement. The large variety of image contrasts and & 
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space sampling methods in MRI causes these artifacts to 
appear differently from scan to scan. This, coupled with a 
large range of motion types that occur in vivo, means that 
there is no single motion correction tool that can be applied 
to every motion problem. Instead, there is a toolbox of 
techniques, where different tools are applicable in particular 
situations. In some cases, powerful tools are already avail- 
able; in other cases, tools still need to be developed. None- 
theless, it is likely that developing dedicated tools for 
specific situations will be what eventually solves the motion 
problem in MRI, rather than a single approach that solves 
all problems. 
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